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In this paper we describe methods for determining the cross-sectional area 
function of the human vocal tract from acoustical measurements made at one 
end. The pressure and volume velocity are assumed to obey Webster's horn 
equation, which is valid for frequencies below 8.5 kHz. Acoustical 
properties below 3.5 kHz do not uniquely specify the area function. This 
paper shows how high frequency information may be incorporated into 
the mathematical model in a manner consistent with a priori information 
about the vocal tract. Some results of application of the methods by computer 
simulation are presented. It is interesting to see from the figures that nine 
numbers (namely, length, four formants, and four residues) specify the 
area function quite well for practical purposes. 

I. INTRODUCTION 

In recent years there has been considerable interest in the modelling 
of speech production in terms of the motion of the articulators. This 
interest has stimulated work on the determination of the shape of the 
human vocal tract as a function of the utterance. For frequencies less 
than 3500 Hz, wave motion in the vocal tract is essentially planar, so 
that the shape is effectively specified by the cross-sectional area as a 
function of distance from one end of the tract (say from the glottis). 

During the past two decades X-ray techniques have been used to 
determine these area functions. These techniques suffer from two major 
drawbacks: (?) In order to keep the exposure to X-rays within safe 
dosage limits, only a small number of measurements can be made on 
any one subject; (it) The interpretation of X-ray data is a complex 
and difficult art, and a number of assumptions must be made in order to 
convert this data to area functions. The accuracy with which area func- 
tions are reconstructed is rather limited. 
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In 1965, Mermelstein and Schroeder suggested the new approach 
of inferring the area functions from acoustic information. 1 Under the 
usual assumptions of lossless plane wave propagation, they showed that 
if the area function A (x) of a vocal tract of length I is of the form 

log A(x) = log A + Z a » cos ~ r W 

n-l * 

then in the limit as a n -» for all n, the nth eigenfrequency (with the 
tract closed at x = and open at re = 1), is given by 

Xn = A „(l - |a 2n _x) (2) 

where X „ is the nth eigenfrequency of the uniform tract (a n = 0, n = 
1, • • •)• Likewise, to the same approximation the nth eigenfrequency 
with the tract closed at both ends, is given by 

/i n = M0n(l - 2fl2n). (3) 

Using equation (2) to obtain a 2n+ i from co n , Mermelstein and Schroeder 
obtained antisymmetric approximations to area functions from a knowl- 
edge of formant frequencies alone. This work was extended by Schroeder 2 
and Mermelstein 8 to include the even-order coefficients [using measured 
values of the poles and zeroes of the input admittance at the lips, which 
correspond respectively to the X's and n's in equations (2) and (3)]. 
An extension was also made by Mermelstein 8 who devised an iterative 
algorithm to compute the first 2m coefficients in equation (1) from a 
knowledge of (Xi • • • X m , mi ■ • • Mm), when the perturbations of the 
eigenvalues are too large for first order perturbation theory to be ac- 
curate. Another iterative scheme has been obtained by J. Heinz, by 
applying perturbation theory to tracts of arbitrary shape. 4 

These methods are applications of very general techniques (namely, 
perturbation theory and steepest descents) which do not make use of 
the special characteristics of the problem at hand. They also leave 
unanswered certain mathematical questions such as the convergence 
of the iterative procedures and uniqueness of the solution. 

In this paper, we describe two (noniterative) methods for computing 
the area function from acoustical data. Apart from clarifying the 
physical and mathematical aspects of the problem, these methods 
provide solutions in a form suitable for analyzing the sensitivity of 
the reconstructed area functions to inaccuracies of the data. They 
also enable us to answer such basic questions as: "What tube has all 
but a finite number of eigenvalues identical to those of a uniform 
tube?" 
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In Section II we introduce the wave equation and Webster's horn 
equation and list the basic properties of the solutions and eigenfunctions 
of the horn equation under homogeneous boundary conditions. 

In Section III we present a method for computing the area function, 
based upon the factorization of the kernel of an integral operator which 
transforms solutions of the horn equation to the solutions of the equa- 
tion for a uniform tract. The existence of such a transform was proved 
by Marchenko, 5 and the transform has been used in the solution of 
the inverse-Sturm Liouville problem by Gelfand and Levitan. 6 

In Section IV we present an alternative method for computing the 
area function based upon the solution of an integral equation whose 
kernel is the driving point response to an impulse at one end of the 
tract. This integral equation was introduced without derivation by 
Krein in a paper on an application of his theory of extensions of positive 
definite kernels. 7,8 Our derivation is physically motivated and uses only 
elementary theory of forced motion of a second order system. 

In Section V we present preliminary results of an application of 
our methods to the determination of vocal tract shapes and a comparison 
with X-ray derived data. Figures 3 through 9 show the results of these 
computations. 

II. MATHEMATICAL PRELIMINARIES 

For a tube of variable cross-sectional area A(x) the equations relat- 
ing acoustical pressure p and volume velocity V are 

(4a) 



§p = 


P dV 


dx 


A(x) dt 


dV _ 
dx 


A(x) dp 
pc 2 dt 



(4b) 

under the assumption of lossless plane wave propagation in the tube. 
These assumptions are accurate for the vocal tract for frequencies 
up to about 4 kHz. For convenience we will choose units such that 
the velocity of sound c = 1, the density of air p = 1 and the length 
of the tube is ir. Then elimination of V in equations (4a) and (4b) 
gives 

s^wg-^g ° = * s * ; (5) 

and for sinusoidal time dependence, such that p = <}>(x, X)e ,x ' ; the 
function <t>(x, X) satisfies 



1198 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1970 



d_ 

dx 



A(p .) a <t>fa x ) + \*A(x)<f>(x, X) = O^x^t (6) 

which is Webster's horn equation. Throughout this paper it will be 
assumed that A(x) > (0 ^ x ^ ir), A(0) = 1, and that A(x) has 
continuous first and second derivatives except at a finite number of 
points in [0, v]. At the points of discontinuity A(x) and its first two 
derivatives are assumed to have finite right and left limits. Under these 
conditions on A (x) the following lemma holds. 

Lemma 1: The solution of equation (6) satisfying the initial conditions 

*(0, X) = 1, *'(0, X) = (7) 



exists, and 



where 



[A(*)]V(*. X) - [A (*)]ty(*, X) | ^ \ (8) 



^ if < °o , £ x £ v, 

and }f/(x, X) is the solution of equation (6) with the same initial conditions 
and A(x) replaced by a canonical shape A (x). The function A (x) is 
such that A o (0) = A(0) = 1, A (x) is constant everywhere in [0, ir] except 
at points of discontinuity of A(x) where it jumps by the same factor as 
does A{x). 

The proof of this lemma is given in Appendix A. 

The solutions of equation (6) satisfying the initial conditions (7) 
become eigenfunctions if they satisfy some homogeneous boundary 
condition at x = ir. These eigenfunctions and eigenvalues have well 
known properties which for the specific case <J)(t, X) = we summarize 
in the following lemma. 

Lemma 2: If A(x) satisfies the conditions described above, then there 
exists a sequence X,- (the eigenvalues) satisfying 

(?') X, > 0, X,- — » oo as i — > » ; 

(ii) 4>{x, X,) are solutions of equation (6) satisfying the initial conditions 
(7) and the condition <f>(ir, X.) = 0; (8a) 

(m) / A(x)<f>(x, \i)<f)(x, X,) dx = 0, iV j, 

= «*, i = h (9) 
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with 

< a; < oo ; 

(iv) 4>(x, X,) are complete in (he space L 2 (0, *■) o/ sguare integrabk 
junctions. 

An immediate consequence of Lemmas 1 and 2 is the following corol- 
lary. 

Corollary 2.1: If n, is the sequence of eigenvalues for the canonical tube 
A (x), with the conditions atx = and x = t as in Lemmas 1 and 2, then 

X,. = Mi . + 00 (10) 

and the a] of Lemma 2 satisfy 

a] = j" A 9 (x)jftz t m) dx + o(|) 



(ID 



- ! °l? 



w/iere ^ is as defined in Lemma 1. 

Finally we will require the following lemma. 
Lemma 3: There exists a function H(x, t) such that 

[A (x)] h +(x, X) - [A(x)Mx, X) + P #(*, 0[A0)ft(l, X) d*. (12) 

This can be proved by substituting equation (12) into equation (6). 
After trivial, but involved, algebra (see Appendix B) it turns out that 
for (12) to be true, H(x, t) must satisfy the following: 

d' Hjx, t) d*H(x. t) , {[4(01*1" „ (r ft = n (n) 

dx 2 d? " + fl^l 3 K } 

^.-o-4r i w i *- ,iAwiiru (u) 



([.4(/)]M'//(.r,/)|, =() -^|^ 



= 0. (15) 



i =.') 



The theory of partial differential equations guarantees the existence 
of a solution to equation (13) under the boundary conditions (14) 
and (15). 
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III. DERIVATION OF A(x) FROM THE SPECTRAL FUNCTION OR TWO SETS 
OF EIGENVALUES 

The spectral function is defined as a staircase function of X with jumps 
of a* at X, (i = 1, • • •)• Thus to say that the spectral function is known 
is equivalent to saying that the pairs (X,- , a]) i = 1, • • • are known. 
Appendix C shows that if \ t (* = 1, • • • ) are the eigenvalues of the same 
tube for the conditions tf>'(0, X) = 0, <f>(0, X) = 1, a<t>(ir, X) + ty'Or, X) = 0, 
(6 ?* 0), then a knowledge of the pairs (X,- , X,-), i = 1, • • • specifies 
the spectral function. Also, in Section IV it will turn out [see equation 
(27), with x = 0] that X, is the ith pole of the driving point impedance, 
and l/2a< the corresponding residue. 

We now derive A (x), given the spectral function. In cases where A (x) 
has continuous first and second derivatives the spectral function 
suffices to uniquely determine A(x). In cases where A(x) has a finite 
number of discontinuities, the locations and magnitudes of the jumps 
are also assumed to be known. Note that equation (12) may be written 
in symbolic form as 

[A (x)H(x, X) = (/ + H)<J>(x, X)[A(x)] J (16) 

where J is the identity on (L 2 [0, x]) and H the integral operator such 
that 

g = Hf** g(x) = f H(x, t)j(t) dt. (17) 

Jo 

Define the operator U which takes a square summable sequence of 
real numbers f it i = 1, 2, • • • to a square integrable function j(x) 
[on (0, it)] defined as 

1(x) = t,1 t [A(x)Mx,\M. (18) 

Define the adjoint operator U* which takes a square integrable function 
f(x) to a square summable sequence /, given by 

/•• = \ f Kx)<p(x, \<)[A(x)] h dx i = 1, ■ • ■ . (19) 

ati Jo 

Let R and R* be defined analogously to U and U*, with [A(x)]*<j)(x, X.) 
replaced by \f/(x, X.-XAo)* in equations (18) and (19). 

Then 

R = (/ + H)U 
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and 

RR* = (I + H)UU*(I + H)*. (20) 

However, from the completeness and orthogonality of the [A(x)]*<f>(x, X,), 
it follows that UU* = I. Thus 

RR* = (I + H)(I + H)*. (21) 

Note that 

[(B8 . _ ^ _ /' /w f (fe^M^M _ «*.y»>) 

•[A (*)Ao(a:)]*<tt (22) 

(because [A (z)]fy(x, m.)/t. is an orthonormal and complete sequence). 

From the asymptotic formulae of Lemma 1, it is seen that the kernel 
of the operator RR* - I is of the Hilbert-Schmidt type [that is, square 
integrable on the square (0 ^ x, t S ■*)>] Therefore, if the X,- , a,- cor- 
respond to those of a tube with appropriate boundary conditions then 
the factorization of equation (21) is always possible. 

This essentially completes our derivation. For the kernel of RR* 
can be constructed if the X.'s and a.'s are known [and in the case of 
discontinuous A(x), the positions and magnitudes of the jumps are 
also known]. The factorization (21) then gives H(x, t). Finally, since 
<l>(x, 0) = 1 and $(x, 0) = 1 (0 ^ x ^ v), equation (12) gives 



[A(x)] J + f H(x, 0[A(0]* At = [AoCaj)] 1 

Jo 



(23) 



which can be solved for [A(x)l J . 

Although, in general, the factorization of equation (18) is difficult, 
we will show in Section V an effective method of computation when 
all but a finite number of X.'s and a,-'s are identical to the corresponding 
Pi's and y/s. 

IV. DERIVATION OF THE AREA FUNCTION FROM THE INPUT IMPEDANCE 

In this section, for simplicity, the area function and its first two 
derivatives will be assumed continuous. Consider the forced pressure 
response y(x, t) in the tube, due to a unit ramp r(t) of volume velocity 
at x = 0. This may be obtained by including a term 8(x)r(t) on the right 
hand side of equation (4b). The resulting equation for y(x, t) is 

f A(x) d -^ - A (x) %-4 . _ s(x)u(t) (a*) 

dx dx at 
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where u(t) is the unit step. Integrating equation (24) over x from 
to a gives 



*«2 



-*»£! 



- f A(x)Z(x, t) dx 1 t ^ 0. (25) 

c = •'0 

In this equation we have put u(l) d 2 y/dt 2 = Z(x, t); it is the transfer 
impedance (in the time domain) from the input end to the point x. 
However (dy/dx) = at x = because of the boundary condition. 
Also, since the velocity of sound has been normalized to unity, for 
t ^ a the region of the tube beyond x = a is undisturbed. Thus for 
t fS a, A (x) dy/dx = for x ^ a. Thus equation (25) becomes 



/; 



A(x)Z(x, t) dx = 1 ^ t ^ a. (26) 

By expansion in terms of (f>(x, A,) it can be verified that 

*<»■« .flfe.a.f +fr*'),""*'' <S (27) 

Of i-1 O!,- 

where the convergence is assumed to be in the sense of distributions, 
and 4>i , a, , A, are as defined in Section II. 
Let j(t) be a function such that 

[ f(t)Z(x, t) dt = 1 x ^ a. (28) 

Then by substitution into equation (26) it follows that 

f /(/) dt = f A(x) dx. (29) 

Jo Jo 

The interesting duality in equations (26), (28), and (29) enables deter- 
mination of A(x) in terms of Z(0, t) rather than Z(x, t). Multiplying 
equation (28) by A(x)Z(x, s), integrating over x and changing order 
of integration on the left-hand side we get 

f f(t) dt [ A(x)Z(x, t)Z(x, s) dx = f A(x)Z(x, s) dx t < a. (30) 

Jo Jo Jo 

For s < a, the right-hand side equals unity by virtue of equation (27). 
On the left-hand side the integration limits on x can be changed to 
(0, v) since Z(x, t) = for x > t. Then substituting for Z(x, t), Z(x, s) 
from equation (27) and using the orthogonality equation (9) we get 



/; 



,/i\ v^ cos X,s cos X,/ j, - , /01 , 

/(0 Z, —2 -dt = 1, s ^ a. (31) 
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[Equation (31) may also be obtained by multiplying equation (28) 
by [A(x)f and using the linear transformation of Section III.] Defining 
f(t) = /(| t |), \t\ ^ a we note that ](t) satisfies 

[" f { t) 2 ^X, S cosX,/ + sinX, S sinX,-/ ^ = ± \ a \g a (32 ) 

since sin X,i is odd. From an elementary trigonometric identity it then 
follows that 

f f(t)Z(0, | t - s |) dt = 1 s ^ a (33) 

and, from equation (29) 

f /(/) o7 = f A(z) dx. (34) 

Jo Jo 

Thus, if Z(0, (which is the driving point impedance function at x = 0) 
is known, or measured, then solution of equation (33) for each a gives 
the area function. [Note that to get A(x) for x ^ a, Z(0, is required 
for t ^ 2a, as expected from physical considerations.]* 

We close this section by noting that although we have discussed 
the method in terms of measurements made at x = 0, where the bound- 
ary condition corresponds to a closed end, trivial modifications are 
needed if measurements are to be made at an open end. In the latter 
case since the pressure vanishes, dv/dx = 0. Therefore, the same method 
is applicable to the horn equation for volume velocity, with a measure- 
ment of driving point admittance (instead of impedance). The driving 
point impedance (admittance) may, of course, be evaluated from 
measurements at any end with an arbitrary, known, termination. 

V. APPLICATION OF THE METHODS TO DETERMINING VOCAL TRACT AREA 
FUNCTIONS 

As noted in the introduction, the one dimensional Webster's horn 

equation is an accurate description of wave propagation in the vocal 

tract, only for frequencies less than about 3.5 kHz. Hence the X, of 

Section II have no physical counterpart whenever they exceed 3.5. 

We therefore start with the \ { and a, (i = 1, • • • n) as measured data, 

and assume that for i > n, the X, and a, are identical with those of some 

canonical tube. In view of the asymptotic formulae of Section II, 

this assumption is reasonable. 

* For another derivation of equation (33) see Appendix D where it is further 
shown that /(a) = [A(a)]" 2 . 
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We know of no a priori method for choosing the canonical shape so 
as to give the best match between the computed area functions and 
those of the actual vocal tract. For simplicity one might assume the 
canonical tube to be uniform. However the experimental area functions 
published by Fant 9 all show a sharp discontinuity at the epiglottis, 
which suggests choosing a canonical tube with such a discontinuity. 
We have tried both these canonical shapes. 

Once a canonical shape has been chosen, a { , X,- , \p(x, X.) and rf/{x, ju,), 
i = 1, • • • n may be computed. Under the assumptions of this section, 
the factorization of equation (21) can then be carried out in the following 
manner. 

We use the vector notation 

kf(z) = M*, X.)/«i . • ' • . W*, M/«- , 

t(x, M<)/Yl , • ' • . t(X, fin)/7n\ 

k?(x) = Mx, X.O/a, ,-■■,*(*, X„)/« n , 

-$(X, Ml)/7l . • * ' . —t(X, Hn)/jn] 

where k^x) and k 2 (x) are rc-dimensional column vectors and the super- 
script T denotes transposition. Then the kernel of equation (22) be- 
comes k\{x)h*{t), and it is easily seen that H(x, t) has the form 

H(x, t) = k?(aj)h(0 x > t (35) 

where h(t) is some vector function of t. Then equation (18) becomes 

k 1 r (.-c)k 2 (/) = kf(x)h(/) + kf(.r) f h(<r)h T (a) dalk^t) (36) 

as long as X,- 5^ /x< ,i = 1, • • • n. (If X, = /*,- for some i, a slight modifica- 
tion is necessary.) However, from the linear independence of the com- 
ponents of ki(aO it follows that 

k 2 (0 = h(0 + f h(*)h T (<r) da\k x (t). (37) 

Equation (37) can be solved for h(t) by the analog circuit shown in 
Fig. 1, or by an equivalent computer simulation. Also, since equation 
(23) now becomes 

[A(x)f + kf(s) f h(/)[A(01* dt = [A (x)f (38) 

•'0 

the analog circuit of Fig. 2 yields [A (a;)]*. 
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*o 



-►h 



hh' 



Fig. 1 — Analog computer circuit for computing h. 

The results of such computations for some area functions published 
by Fant are shown in Figs. 3 through 9. 

We close this section by noting that if the experimental data is in 
the form of a driving point impulse response, then the simplest procedure 
is to use the method of Section IV [that is, to solve equation (33) for 
various values of a]. We have not computed area functions by this 
method so far, but propose to do so, using impedance tube or other 
experimental data. The limitations due to the inapplicability of the 
horn equation at high frequencies apply to this method as well. The 
effect of low-pass filtering the driving point response is being inves- 
tigated. 

VI. CONCLUSIONS AND DISCUSSION 

The comparison between measured and computed area functions 
of Figs. 3 through 9, indicates that knowledge of the first few (X, a) 
pairs is sufficient to get reasonable estimates of the area function. 
The X's may be obtained directly from the speech output, since they 
can be computed with reasonable accuracy from the formant frequencies. 
The a's on the other hand cannot be computed directly from the speech 
waveform, and impedance tube or other equivalent measurement would 
appear to be necessary. However, the vocal tract has physical con- 
straints which might be reflected in a functional dependence of the 
a's on the X's. The possibility of such functional dependence is being 



(Ao)* ±0 



■(A)' 



KS> <| & 



Fig. 2 — Analog computer circuit for computing A. 
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Fig. 3 — Area functions reconstructed from the first three poles and residues of 
input impedance by the method of Section V using a uniform canonical tube. 
Dashed curves are the X-ray measurements. 

investigated. The sensitivity of the computed area functions to changes 
in the a's is also being investigated. 

If one is willing to make acoustical measurements at the lips, then 
the method of Section IV is the most direct way of computing the area 
function. It has the added advantage that the length of the vocal tract 
need not be known. Some preliminary results on the effect of band- 
limiting the impulse response have been obtained and will be reported 
in a later paper. 

APPENDIX A 

Proof of Lemma 1 

Under the assumptions of this lemma, equation (6) may be trans- 
formed to: 

[[A(x)Mx, y))" + X'[1W1W. X) = {[AixtfVvix, A) (39) 
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except at the points .r, , • • • x k , where A{x) is discontinuous. With 
x = and x k+t = ir, equation (39) gives 

[A(x)Mx,\) 

= Ux, x) + I — \ — ~[M(ft [M0Mt ' X) dl ' 

x { < x ^ x i+x , i = 0, 1, ••• ,k (40) 



vhere 



/,(x, X) = a,(X) cos \x + 6i(X) sin Xz. 



The coefficients a,(X), 6,(X) are to be determined so as to make <p(x, X) 
and A(x)<p'(x, X) everywhere continuous. (The conditions at x - 
give a (X) = 1, 6 (X) = 0.) Clearly for every X there exists a bound 
m,(X) = sup | [A(x)] h <p(x, X)|, x,- ^ x ^ z i+1 . Then from equation (40), 




Fig. 4 — Same as in Fig. 3, except the canonical tube was chosen with a dis- 
continuity at the epiglottis. 
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m t 



[Am 



ii/f 



[MW 



dt 



i<(X) ^ I /,(*, X) | + £ m..(X) jf | sin \(x - I) \ 

^ [a%\) + tf(X)]* + Mm,(X)/X (41) 

where M is a bound on the integral for all i. Thus for X > 2M, 
| [AfeOlW, X) | ^ 2[o?(X) + 6?(X)]* 
= 2y,(X). 



(42) 
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Figs. 5-9 — Area functions reconstructed from the first four poles and residues: 
(a) the reconstruction using a uniform canonical tube, and (b) the reconstruction 
with a discontinuous canonical tube as in Section V. Dashed curves are X-ray 
measurements. 
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Then from equation (40), 

Xi < x < Xt+i (43) 



[A(x)Mx, X) = /,(.r, X) + 7,(X) x C,Ct) , 



where c,-(x) is bounded. Differentiating equation (40) and using a 
similar argument gives, 



[A (a;)]* d<p(x, X) _ 1 d/,-(.r, X) 



X dx X dx 

with d,(x) bounded. Denning 

"[A(.r,.+)/A(.T,.-)] 



- 7,(X) di(x)/\, x { < x < x i+J (44) 



k, 







[A(.r 1 -)/A(.r..+)] i . 



and 



r. 



R, = 



cos Xx, sm \x { 
Lsin \Xi — cos \x, J 
the continuity conditions at x t give: 



o,-(X) 
L6,(X)J 



= R.k.R, 

+ 



o,-,(X) 

L&,- 1 (x)J 

7«-i(X) 



R,k, 



ldi-j(Xi) 



7,00 



R, 



C(x-,.) 

Mxdj 



(45) 



Since the norm of R, is unity and of k,- finite it follows upon taking 
the lengths of the vectors on either side of equation (45) that for large 
enough X, 7i(X) ^ K'yi-i(\), for some finite constant K' . Since 7o(X) = 1, 
it follows that 7,(X) is bounded for all i, as X — » °o . Then from equation 
(45) 



Lbi_ 



= Rik.R; 






+ 



(46) 



However, if a, , b t satisfy equation (45), then/,(.c, X) = [A ,(x)\ k \f/(x, X) -f- 
0(1/X) for Xi ^ x ^ x, + i , therefore from equation (43) 



[A(x)Mx) = [A Q (x)\H(x) + 



X 



(47) 



with c(x) bounded. 
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APPENDIX B 

Existence of the Linear Transformation of Equation (12) 

In this appendix we prove the existence of the linear transformation 
of equation (12). Consider first the region ^ x ^ x x , where Xi is 
the first point of discontinuity of A(x). Then with y(x) = [A(x)f4>(x, X) 
and q(x) = [[A(,xyf)"/[A(xyf s equation (6) becomes 



y"(x) = -\'y(x) + q(x)y(x). 
Consider the function 



3l(x) = f H(x, t)y(t) dt. 



(48) 



(49) 



Then 



9l"(aD = jr [H(x, x)y(x)] 



+ *£** 



+ / 



1 d 2 H(x, t) 
dx 2 



y{t) dt. (50) 



After integrating twice by parts we arrive at the identity 



£ 



d 2 H(x, t) , A ,. dU(x, t) ,, 



-*(*,»§ 



+ f H(x, t) ^ dt. (51) 



Substituting for d 2 y/dt 2 from equation (48) into equation (51) and 
adding equation (50) we get 



ttww;[^-^+'«M 



■y(t)dt + 2y(x)j-H(x,x) 



-y(t) 



dH(x, t) 
dt 



dy(t) 



+ H(x,t) 

1-0 «*' 



(52) 



If now H(x, I) satisfies the differential equation (13) with boundary 
conditions (14) and (15), then equation (52) shows that 31(a) + y(x) 
is some linear combination of cos \x and sin Xx. Matching of boundary 
conditions at x = then shows that 

Sl(a;) + y(x) - cos \x = [A (x)]ty(x, A), ^ x ^ x, . (53) 
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The proof may be extended to x > z, by a similar procedure. Thus, 
for example, in the range x x ^ x ^ x 2 , H(x, t) must satisfy the dif- 
ferential equation (13), the boundary conditions (14) and (15) for 
x > Xx , and boundary conditions at x = x x imposed by the continuity 
requirements on </»(z, X), f(x, X), A(x)<t>'(x, X) and A (x)4>'(x, X). 

APPENDIX C 

Spectral Function from Two Sets of Eigenvalues 

We outline a method of getting a spectral function from two sets of 
eigenvalues. Let <p{x, X) be the solution of equation (6) such that 

ipfa X) = a, A(*W(t, X) = j8 (54) 

for every X. Let \l/(x, X) be the solution such that 

$( V , X) = y, A(tW(t, X) = 5. (55) 

Let X? , X 2 , , • • • be the values of X 2 for which a<p(0, X) + &A(0V(0, X) 
= and let jx? , fh , "' be the values of \i for which a^(0, n) + 
bA(0)f(0,iO = O.Let 

m(X) = wKO, X) + 6A(0)/(0, X) (5b) 

Then the zeroes of m(X) are fi x , p» , • • • and the poles are Xi , X 2 , • • • . 
If X(x, X) is any solution of equation (6), then it is easily shown that 

(X 2 - X 2 ) P A(x)X(x, \)<p(x, X n ) dx 

Jo 

= A(x)[X(x, \)<p'(x, K) - <p(x, K)X'(x. X)] * . (57) 

Choosing X(x, X) = fix, X) - m(\)<p(x, X) in equation (57) and using 
the boundary conditions on ${x, X) and <p(x, X), we get 

(X 2 _ X 2 ) P A(x)[+(x, X) - m(\)<p(x, X)] dx = fa - aS (58) 
Jo 

for all X. As X — > X„ 

<xl = [' A(x)v\x, *») dx = lim (aS - fa)/[(\ 2 ~ X 2 )m(X)]. (59) 

Jo x-x» 

Thus, given X, , X 2 , and hi , n 2 , one obtains m(X) and hence a t , a a , • • « , 
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APPENDIX D 



Derivation of Integral Equation (33) 

We give here a derivation of the integral equation (33) based upon 
the results of Section III. For simplicity we will assume that A(x) 
[and hence A (x)] has no discontinuities. Then from equation (23) 

(7 -f- H)A*(x) = u(x) (60) 

where u(x) is equal to 1 for all x > 0. Thus if j(x) is a function such that 

[(/ + Hy l u]{x) = g(x) (61) 



then 



Notice that if 



then 



r 



r g 2 (x) dx = f A(x)dx. (62) 

Jo Jo 

(I + H){I + H)* = I + K (63) 

A(x) dx = f [(/ + Ky l u](x) dx 

/(.t) dx (64) 



where j{x) satisfies the equation 

[(/ + K)j\(x) = u(x). (65) 

The kernel of I + K is recognized as that of equation (31) with a = ir. 
Equation (33) therefore follows (for a = t) from the symmetrization 
of j(x), exactly as in Section IV. However, the argument given here 
is independent of the length ir, which may be replaced by a. 
Using equations (60), (63) and (65), we have 

(/ + #)*/(.) = [A(-))K (66) 

Therefore from equation (17) 

/(a) = [A(a)]K (67) 
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